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Abstract — This  paper  covers  the  second  part  of  the  analysis  of 
data  recorded  by  the  surface  wave  (SW)  over-the-horizon 
(OTH)  WEllen  RAdar  (WERA).  Data  were  collected  by  two 
WERA  systems,  on  May  13th  2008,  during  the  NURC 
experiment  in  the  Bay  of  Brest,  France.  The  principal  aim  of  this 
work  is  to  provide  an  accurate  characterization  of  the  spectral 
components  of  the  received  signal.  Secondly,  this  information  is 
exploited  in  order  to  provide  a  simple  and  reliable  spectral 
modeling  tool.  For  this  reason,  auto-regressive  (AR)  models,  also 
known  as  linear  prediction  (LP)  models  have  been  investigated. 
Our  results  show  that  at  long  distances,  when  the  clutter-to- 
noise  power  ratio  (CNR)  is  small,  the  main  components  of  the 
spectrum  can  be  reasonably  described  by  an  AR(12)  model,  with 
a  good  compromise  between  accuracy  and  simplicity.  As  the 
CNR  increases  higher-orders  are  instead  to  be  preferred. 

I.  Introduction 

High-frequency  surface- wave  radar  (HFSWR)  systems  are 
particularly  attractive  for  remote  sea-state  sensing  by  virtue  of 
their  over-the-horizon  scanning  capability  [1].  Low-power 
HFSWRs  can  be  found  operating  from  the  coast  or  on  board 
of  ships.  This  is  the  case  of  the  WEllen  RAdar  (WERA), 
developed  at  the  University  of  Hamburg  and  mainly  devoted 
to  remote  sensing  purposes  [2],  [3].  By  exploiting  the  good 
conductivity  of  sea  water  in  the  lower  HF  band  (3-15  MHz), 
HFSWRs  are  successfully  used  for  detecting  vessels  and  ships 
at  long  distances  in  the  context  of  maritime  surveillance  and 
borders  monitoring.  In  fact,  they  are  able  to  solve  the  problem 
of  line-of-sight  (LOS)  coverage  typical  of  microwave  radars 
and  do  not  suffer  from  the  periodic  coverage  of  satellite  based 
sensors  [4],  [5].  However  high-power  systems  are  quite 
complex  and  expensive,  so  even  growing  interest  is  focused 
on  low-power  HFSWRs,  as  WERA. 

The  contribution  of  HFSW  sea  clutter  is  produced  by 
specific  spectral  components  of  the  surface-height  wavefield. 
First-order  Bragg  scattering  is  due  to  those  ocean  waves  of 
half  of  the  radar  wavelength  which  travel  towards  and  away 


from  the  radar  site.  The  Doppler  spectrum  of  the  backscattered 
signal  then  contains  two  lines,  corresponding  to  the  phase 
velocities  of  the  scattered  ocean  waves.  These  frequencies 
often  deviate  from  the  theoretically  known  values  in  non¬ 
moving  waters.  This  phenomenon  is  due  to  an  underlying 
surface  current  [6],  [7].  Furthermore,  in  addition  to  the  Bragg 
scattering  returns,  a  more  complicated  feature  is  referred  to  as 
second-order  scattering. 

The  working  range  of  a  HFSW  radar  depends  not  only  on 
the  attenuation  of  the  electromagnetic  wave  during  the 
transmitter-target-receiver  path,  but  also  on  the  scattering 
strength  of  the  target,  atmospheric  noise  and  noise  due  to  radio 
interference  [1].  In  fact,  the  ionosphere  is  prone  to  propagate 
interference  at  long  distances,  especially  at  night. 
Furthermore,  external  interference  from  natural  and  man-made 
sources  typically  masks  the  entire  range-Doppler  search  space 
and  is  characterized  by  a  spatial  covariance  matrix  that  is 
time-varying  or  non- stationary  over  the  coherent  processing 
interval  (CPI).  This  physical  phenomenon  may  arise  from  a 
number  of  causes.  For  instance,  the  dynamic  properties  of  the 
ionosphere  layers  propagating  the  HF  interference,  the 
variation  in  geometry  between  radar  receiver  and  interference 
sources  and  the  impulsive  nature  of  the  sources  [1],  [8].  For  all 
these  reasons,  OTH  systems  still  represent  a  great  challenge 
for  the  scientific  community  [2]. 

Regarding  ship  detection,  the  task  is  to  resolve  targets  in 
the  temporal  or  in  the  Doppler  frequency  domain  from  the 
same  background  clutter  exploited  for  sea-state  sensing.  The 
two  problems  are  thus  in  many  ways  complementary.  In  fact, 
the  presence  of  clutter  is  unwelcome  as  far  as  ship  detection  is 
concerned,  while  the  presence  of  ship  returns  can  negatively 
affect  the  extraction  of  oceanographic  parameters.  For  all 
these  reasons,  in  the  past  years  there  has  been  much  interest  in 
trying  to  develop  new  spectral  techniques  for  modeling  the 
return  from  the  sea,  with  the  ultimate  goals  of  enhancing 
performances  in  both  target  detection  via  clutter-suppression 
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techniques  [9]  and  sea  current  sensing  [10],  by  means  of 
knowledge-based  receiver  configurations. 

This  paper  deals  with  the  spectral  analysis  of  sea  clutter 
using  signals  received  by  two  concurrently  operating  WERA 
systems  and  it  is  organized  as  follows.  A  brief  overview  about 
WERA  and  the  NURC  experiment  in  the  Bay  of  Brest 
(France)  is  provided  in  section  II.  In  Section  III  a  qualitative 
study  about  the  main  spectral  components  of  the  HE  signal  is 
carried  out,  while  in  Section  IV  the  AR  modeling  of  sea  clutter 
spectrum  is  presented.  Conclusions  are  instead  drawn  in 
Section  V,  in  addition  with  some  guidelines  about  future 
work. 

II.  WERA  System  and  Experiment  Description 

The  WERA  system  is  essentially  a  shore  based  remote 
sensing  system  for  monitoring  ocean  surface  currents,  wave 
spectra  and  wind  direction.  Its  development  has  been  carried 
out  in  1995  at  the  University  of  Hamburg  by  K.  W.  Gurgel 
[2],  [3].  Data  were  collected  during  the  NURC  experiment  at 
Brest  (Brittany,  France),  on  the  13th  of  May  2008  [11].  The 
two  WERA  systems,  respectively  located  at  Brezellec 
(Latitude  48°  4’  8”,  Longitude  4°  40’  0”)  and  Garchine 
(Latitude  48°  30’  10”,  Longitude  4°  46’  32”)  are  owned  by 
the  French  Service  Hydrographique  et  Oceanographique  de  la 
Marine  (SHOM)  and  operated  by  the  company  SAS  ActiMar. 
They  were  located  parallel  to  the  coast  and  aligned  towards 
the  true  North  with  angles  of  20°  and  339°,  respectively.  Data 
records  occurred  every  twenty  minutes  at  each  receiver  site 
(i.e.  Brezellec  and  Garchine)  and  last  for  8  minutes  and  52 
seconds  (corresponding  to  a  CPI  of  532  seconds),  for  a  total  of 
2048  measurements  with  a  chirp  duration  of  0.26  s,  suitable 
for  sea  surface  variability.  Both  the  WERA  systems  worked  at 
frequencies  comprised  between  12.190  and  12.565  MHz  [12]. 
Bandwidth  was  set  to  100  kHz  (i.e.  the  range  resolution  is  1.5 
km).  For  each  chirp,  100  range  samples  were  collected  (for  a 
maximum  range  of  1 5  0  km)  [13]. 

III.  Spectrae  Anaeysis 

This  paper  is  devoted  to  analyze  and  model  the  spectral 
features  of  the  received  HE  signal.  Since  the  receiver  is  a  16 
element  linear  array  of  quarter-wavelength  monopoles, 
initially  a  qualitative  analysis  was  carried  out  for  each  antenna 
separately  [14].  Before  the  statistical  analysis,  data  have  been 
beamformed  by  means  of  a  simple  rectangular  window  [15]. 
Two  different  data  records  have  been  considered,  as 
representative  of  all  the  files:  Brezellec  dataset  no.  8  (since 
02:30:00  through  02:38:52  UTC)  and  Garchine  dataset  no.  6 
(since  01:40:00  through  01:48:52  UTC).  The  aim  is  to  focus 
our  attention  on  the  most  important  features  of  HE  signals. 

A.  Sea  clutter  analysis 

A  two-dimensional  color  plot  describes  the  behavior  of  the 
estimated  power  spectral  density  (PSD)  in  [dB]  as  a  function 
of  the  range  cells  (rows)  and  of  the  Doppler  frequency 
(columns),  as  shown  in  Fig.  1  for  the  8^^  file  recorded  at 
Brezellec.  A  512  samples  window  (corresponding  to  133.12  s) 
and  an  azimuth  angle  of  15°  have  been  considered.  For  the 
majority  of  range  cells,  two  lines,  which  are  related  to  the 
advancing  (positive  frequency)  and  receding  (negative 


frequency)  waves  are  easily  recognized.  The  frequency 
separation  A/’  between  the  two  first  order  Bragg  peaks  can  be 
evaluated  using  the  following  formula  [1]: 

A/  =  2  (1) 

V  7t'C 

where  g  is  the  acceleration  due  to  gravity,  c  is  the  speed 
of  light  and  is  the  centre-band  frequency.  The  Doppler 
separation  corresponds  to  about  0.72  Hz  for  all  the  operative 
frequencies  comprised  between  12.190  and  12.565  MHz. 
Second-order  peaks  are  well-defined  in  the  first  forty  range 
cells.  In  this  interval  sea  clutter  level  dominates  both  targets, 
and  noise.  Anyway,  as  we  move  away  from  the  receiving 
radar,  the  CNR  decreases  and  the  two  Bragg  peaks  are  less 
defined,  with  the  final  result  that  background  noise  overcomes 
them.  In  this  region  sea  clutter  tends  to  become  a  white 
random  process. 


Figure  1.  Brezellec  dataset  no.  8:  clutter  PSD  in  [dB]  for  an  azimuth  angle 

of  15°. 

Together  with  sea  clutter,  especially  for  high  CNR  levels, 
the  contribution  of  land-scattered  echoes  around  the  dc 
frequencies  can  be  significant.  In  addition,  the  two  Bragg 
peaks  are  subject  to  a  frequency  modulation.  Rapid  changes 
(e.g.  turbulences)  in  the  sea  waves  movement  around  rock 
cliffs  lead  to  an  evident  broadening  of  the  two  Bragg 
frequencies  [14],  as  shown  in  Fig.  2  for  the  6^^  file  recorded  at 
Garchine.  This  phenomenon  happens  for  all  the  files  collected 
at  Garchine,  thanks  to  the  presence  of  large  islands  (e.g.  He 
d’Ouessant)  at  relatively  small  distance  from  WERA. 

B.  Interferences  analysis 

In  HFSWRs,  sea  clutter  is  not  the  only  type  of  signal  we 
can  observe,  since  systems  are  prone  to  a  number  of 
interference  sources,  both  natural  and  man-made.  The  former 
usually  consist  of  large  returns  (horizontal  lines)  that  cover  a 
large  portion  of  the  Doppler  space.  In  this  specific  case  they 
extend  from  about  0.20  to  1.60  Hz  and  manifest 
approximately  around  range  cell  85  (about  126  km  away  from 
the  coast),  as  shown  in  Fig.  1.  These  interferences  are 
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responsible  of  the  heavy-tailed  behaviour  of  clutter  amplitude 
described  in  [11].  They  are  principally  due  to  unwanted 
propagation  modes  through  the  ionosphere  and/or  meteor 
trails  echoes. 

The  second  type  of  interference  is  represented  by  radio 
frequency  interference  (RFI),  as  shown  in  Fig.  2.  These 
returns  mask  both  sea  clutter  and  ship  echoes  in  even  larger 
portions  of  the  range-Doppler  (RD)  space.  In  fact  they  appear 
in  certain  Doppler  frequencies  intervals  as  vertical  lines 
extending  for  all  the  defined  range.  For  overcoming  this 
serious  problem  an  algorithm  was  proposed  in  [12]  for 
estimating  and  then  removing  RFI  noise. 
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Figure  2.  Garchine  dataset  no.  6:  elutter  PSD  in  [dB]  for  an  azimuth  angle 

of  10°. 


IV.  Auto-Regressive  Modeeing 
A  random  process  X[n)  is  an  auto-regressive  (AR)  model 
if  it  can  be  described  by  the  following  expression: 


X{n)  =  ^a,,X{n-k)  +  W{n), 


(2) 


where  (for  k  =  are  the  model  parameters 

and  W[n)  is  a  white  Gaussian  process  with  zero  mean  and 
variance  .  The  PSD  (/)  process  is  given  by: 


5r(/)= 


l-^a,exp{-72M/r} 


-,for|/|<l/2r.  (3) 


kp  dp  2 

^  T 

ap  p  j  , 

(5) 

II 

'■uur, 

(6) 

(0)  (1) 

-  r,{P-2) 

o 

-  ryp-\) 

,  (7) 

■yp-\)  ryp-2) 

o' 

depend  onr^  (m)  evaluated  from  the  set  of  data: 

rx{m)  =  E^X (n-m)X*  .  (8) 

Then,  the  estimated  variance  of  the  driving  process  noise  is 
obtained  as  follows: 


k^l 

Finally,  we  substitute  and  ,  obtained  respectively 
with  (4)  and  (9),  in  (3). 

A.  Modelling  accuracy 

Let  us  consider  dataset  no.  8  recorded  at  Brezellec.  Some 
preliminary  results  are  shown  in  figures  3  and  4,  respectively 
for  range  cells  2 1  and  87  and  for  a  common  azimuth  angle  of 
0°.  They  show  the  comparison  between  the  estimated 
periodogram  (blue  line)  and  the  PSD  obtained  by  four  AR 
models:  AR(6),  AR(12),  AR(24)  and  AR(48),  represented  by 
the  yellow,  green,  magenta  and  red  curves  respectively. 

The  main  result  is  that  for  near  range  cells,  for  which  the 
CNR  is  considerable,  second-order  peaks  can  be  optimally 
described  only  by  higher-order  models,  such  as  the  AR(48) 
which  shows  to  be  the  most  effective  in  terms  of  modeling 
accuracy  of  second-order  scattering  returns  (Fig.  3).  Low- 
order  AR  models  are  able  to  follow  only  the  main  features  of 
the  signal,  but  not  second-order  contributes. 

Anyway,  for  distant  range  cells,  second-order  Bragg 
scattering  becomes  more  negligible,  thus  allowing  us  to  adopt 
a  simpler  AR(12)  model,  which  guarantees  a  good 
compromise  between  accuracy  and  simplicity  (Fig.  4). 
Conversely,  low-order  AR  models  tend  to  generate  broad¬ 
band  first-order  Bragg  contributes,  as  done  by  the  AR(6),  and 
even  in  this  case  are  not  suited  at  all. 


Chosen  a  desired  order  P,  the  coefficients  apj^  are 

estimated  by  solving  the  system  of  Yule- Walker  set  of 
equations  expressed  by  [15],  [16]: 

V  =  ’  (4) 


where: 
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Figure  3.  Brezellec  dataset  no.  8:  elutter  PSD  vs  AR  modelling  in  [dB]  for 
range  e ell  21  and  azimuth  0°. 


Figure  4.  Brezellee  dataset  no.  8:  elutter  PSD  vs  AR  modelling  in  [dB]  for 
range  eell  87  and  azimuth  0°. 


B.  Choice  of  the  order  P 

A  common  procedure  for  estimating  the  most  suitable 
order  consists  in  evaluating  the  root  mean  square  error 
(RMSE)  of  the  estimated  PSD: 


P-.RMSE{^P)  =  mm[RMSE{p)].  (11) 

This  procedure  is  carried  out  for  all  the  cells  in  RA  space 
(i.e.  =100  and  =121,  where  the  numer  of 

azimuth  cells  by  which  we  divided  the  angle  of  view  after 

beamforming).  The  behaviour  of  P  ,  based  on  the  test  in  (1 1), 
is  plotted  against  all  RA  cells,  as  shown  in  Fig.  5.  The 
estimated  order  is  generally  smaller  than  10.  However,  in 

some  areas  P  is  high,  as  in  the  high-rigth  side  of  the  figure  5 
between  40°  and  60°.  Here  clutter  is  characterized  by  several 
peaks  in  the  Doppler  frequencies  spectrum,  thus  meaning  that 
the  optimal  AR  order  is  reasonably  high. 


5  10  15  20  25  30  35  40  45 


Azimuth  ['*] 


Figure  5.  Brezellec  dataset  no.  8:  estimated  P  order  at  the  varying  of  RA 

cells. 

A  further  analysis  is  carried  out  after  averaging  the  RA 
profile  of  P  over  azimuth  cells,  as  shown  in  Fig.  6.  The  best 
choice  lies  between  10  and  15  for  the  majority  of  the  cells 
thus  confirming  the  hypothesis  that  the  AR(12)  model  is  a 
reasonable  good  choice.  Conversely,  higher-order  models 
should  be  used  for  modeling  near  RCs,  while  lower-order 
models  for  distant  cells.  The  peak  around  RC  80  (about  120 
km  away  from  WFRA)  is  caused  by  the  presence  of  a  strong 
echo  near  zero  Doppler  frequencies,  as  shown  in  Fig.  1.  This 
peak  appears  for  the  majority  of  azimuth  cells. 


1/2  2 

RMSE{p)^\\[s,{f)-ST{f)\df,  (10) 

V-1/2 

where  is  the  empirical  PSD  of  the  process 

estimated  from  the  data  (the  periodogram)  and  (/)  is  the 
AR  power  spectral  density,  for  a  given  set  of  orders  (i.e. 

=  2,4,6,...,48).  The  choice  of  P  ,  for  a  given  RA  cell,  is 
done  considering  the  value  of  p  which  minimizes  the  RMSF. 
The  decision-test  is  given  by: 


Concluding,  we  can  state  that  the  AR(12)  is  a  good  model 
both  in  terms  of  simplicity  and  accuracy.  Figures  7  and  8 
respectively  show  the  range-Doppler  (RD)  behavior  in  [dB]  of 
the  periodogram  and  the  PSD  estimated  by  means  of  the 
AR(12)  model,  for  an  azimuth  angle  of  0°.  As  we  can  observe, 
the  main  features  of  the  signals  are  almost  always  preserved. 
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Figure  6.  Brezellec  dataset  no.  8:  estimated  P  order  as  a  funetion  of  range 
[km],  data  have  beeen  averaged  over  azimuth. 
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Figure  7.  Brezellee  dataset  no.  8:  empirie  PSD  in  [dB]  for  an  azimuth  angle 

of0°. 
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Figure  8.  Brezellee  dataset  no.  8:  PSD  in  [dB]  modelled  by  an  AR(12)  for 
an  azimuth  angle  of  0°. 


C.  Analysis  of  the  poles 

In  this  Section,  the  poles  of  the  AR  model  are  analyzed 
and  their  distribution  in  the  complex  plane  is  discussed  [9]. 
Two  cases  are  presented,  both  referring  to  dataset  no.  8 
recorded  at  Brezellec:  the  AR(6)  and  the  AR(IO),  as  shown  in 
figures  9  and  10,  respectively.  In  both  cases  the  poles  of  the 
AR  models  (depicted  by  blue  “x”)  are  considered  for  all  the 
RCs  between  5  and  100.  For  completeness  the  two  phases 
corresponding  to  the  advancing  and  receding  Bragg  waves  are 
reported  as  well  (red  dashed  lines). 

As  expected  all  the  poles  are  inside  the  unit  circle,  thus 
confirming  the  stability  of  the  model.  Two  main  accumulation 
clouds  appear  very  close  to  the  unit  circle  [9].  This  means  that 
the  HF  signal  has  two  principal  narrow-band  components. 

Usually  this  pair  of  poles  starts  moving  from  the  unit  circle 
towards  the  centre,  where  the  other  poles  accumulate.  This 
phenomenon  is  justified  by  the  fact  that  with  increasing 
distance,  the  contribution  of  sea  back-scattering  becomes  more 
negligible,  as  shown  in  figures  1  and  2.  Therefore,  the  overall 
signal  tends  to  become  a  white  noise  process,  as  confirmed 
also  by  the  estimate  of  the  AR  order  described  in  Fig.  6,  for 
which  the  order  is  generally  smaller  than  5. 

The  remaining  poles  spread  inside  the  unit  circle  with  the 
same  phase  spacing  and  the  same  radius,  except  for  spurious 
contributes  due  to  interference  sources.  The  task  of  these  poles 
consists  in  describing  the  “floor”  part  of  the  spectrum.  With 
increasing  order,  they  move  towards  the  circle,  because  the 
new  spectral  contributes  become  even  more  negligible  if 
compared  to  the  two  Bragg  components  (Fig.  10).  In  addition, 
when  the  order  is  sufficiently  high,  the  almost-dc  component 
of  the  signal  is  described  by  a  pole  on  the  real  axis,  usually 
close  to  the  unit  circle,  especially  for  near  range  cells. 
Concluding,  in  the  case  of  interferences,  new  poles  are  moved 
towards  the  disc,  for  granting  and  overall  good  spectral 
approximation. 
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Figure  9.  Brezellec  dataset  no.  8:  zeros  and  poles  of  the  AR(6)  model,  range 
cells  5-100  and  for  an  azimuth  angle  of  0°. 
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Figure  10.  Brezellec  dataset  no.  8:  zeros  and  poles  of  the  AR(IO)  model, 
range  eells  5-100  and  for  an  azimuth  angle  of  0°. 

V.  Conclusions 

In  this  paper  an  accurate  analysis  of  the  spectral  features  of 
sea  clutter  in  the  HF  band  has  been  carried  out.  Initially  we 
provided  a  qualitative  investigation  of  the  main  spectral 
components.  Several  interesting  features  were  observed.  The 
two  Bragg  lines  which  characterize  the  long-term  sea  clutter 
variations  are  well  defined  for  the  majority  of  cells.  On  the 
contrary,  the  second  order  peaks,  which  are  instead 
responsible  of  the  rapid  changes  in  the  wave  behavior,  are  well 
defined  only  for  near  range  cells.  The  contribution  of  land 
scattering  can  be  observed  as  well  (but  only  for  near  mixed 
sea/land  cells)  in  addition  with  unwanted  signals,  both  man¬ 
made  and  natural.  Variations  in  the  RD  profile  strongly 
depend  on  the  angle  of  view,  after  beamforming.  In  the  second 
part  we  presented  a  common  AR  modeling  technique.  Our 
results  have  evidenced  that  only  for  low  CNR  levels  it  is 
possible  to  describe  the  clutter  signal  with  a  reasonably  low- 
order  AR  model.  For  reducing  the  spectral  parameters  under 
estimations,  we  tried  also  a  different  parametric  model 
composed  of  simple  functions  as  the  Gaussian  and  the 
exponential.  Unfortunately,  the  results  were  not  promising  and 
we  decided  to  omit  them.  The  ARMA  modelling  is  currently 
under  investigation  and  will  be  presented  in  future  works. 
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